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We discuss the properties of a class of saddle point solutions in SU(2) in three dimensions (SU(2) 3 ), exhibiting 
localized peaks in the action. These configurations are generated by deterministic cooling and extremizing algo- 
rithms from analytic configurations. They share some characteristics with cooled and extremized Monte Carlo 
generated lattices. We have investigated physical behavior such as the string tension by averaging over this class 
of saddle point configurations. We have also measured the eigenvalues for harmonic fluctuations around these 
configurations. 



1. Introduction 

Confinement in non-Abelian gauge theories has 
been a very challenging problem to understand 
analytically. Even though there is considerable 
evidence for confinement at weak coupling from 
numerical simulations, the underlying property of 
gauge field configurations giving confinement has 
not been understood. This is largely due to the 
fact that confinement is a non-perturbative phe- 
nomenon, for which we have limited tools for in- 
vestigation. 

Cooling [1,2] (deterministic reduction of the 
action) or extremizing [3] (deterministic reduc- 
tion of the Hessian, S = *£i, x {d S / d A i {x)) 2 ), 
when applied to Monte Carlo generated SU(2) 
lattices, tends to produce configurations with iso- 
lated peaks in the action density. And, remark- 
ably, the string tension sustains its value after 
cooling has decreased the action to less than 1 
percent of its initial value [1]. Since there are 
no stable nonzero solutions in SU(2)3, the quasi- 
stability under cooling is likely due to being close 
to a saddle point solution with a few unstable 
modes. 

Considering the fact that the saddle point can 
be defined precisely, in contrast to a "moderately 
cooled lattice" , extremized lattices serve as a bet- 
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ter staring point for quantitative investigation. 
Using these saddles, we could even try to get a 
valid approximation of the partition function by 
expanding around these saddle points: 

Z = jdu ie - s ^ 

= y I du ie - s ^ 

i 

where JJ^ dUi covers the gauge link configurations 
near a particular semiclassical solution, and is 
the fluctuation determinant near that solution, 
which is given by the product of the eigenvalues 
of harmonic fluctuations. (In the Georgi-Glashow 
model [4], the harmonic fluctuations were decou- 
pled from the particular saddle, = constant.) 

To achieve this, we need to (i) develop a sys- 
tematic way to cover all the saddle points and (ii) 
get the eigenvalue spectrum around each saddle 
point. In this work, we use a starting configu- 
ration to make saddles which shares one char- 
acteristic common to many saddles found from 
Monte Carlo generated lattices. Section 2 de- 
scribes the starting configuration and how we gen- 
erate the saddle point solutions with this config- 
uration. Section 3 discusses the result of string 
tension measurements of these saddle points. Sec- 
tion 4 discusses the eigenvalues of harmonic fluc- 
tuations around these configurations. 
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Figure 1. The starting configuration (Z(2) vor- 
tex) we use to produce saddle-point configura- 
tions. The monopole and anti-monopole (M and 
M) are defined relative to an arbitrary U(l) sub- 
group of SU(2). Wilson loops centered on the 
point P have a value of —1. Two such loops are 
shown by the dotted lines. 



Figure 2. The action density plot of saddle-point 
configurations generated from a Z(2) vortex 



2. Generation of saddle gas 

On cooled lattices, we observed that Wilson 
loops centered on an action peak have trace very 
close to —1. This is the indication of Z(2) flux 
associated with the peak. So we have come up 
with the starting configuration in Fig. 1, where 
Wilson loops around P have trace = — 1. 

Fig 1 shows a monopole-antimonopole pair 
with respect to an arbitrary U(l) subgroup of 
SU(2). Flux from the monopoles is confined to 
the closed tube. Because this configuration is 
compact in size, we can distribute more than 
one of these vortices on a lattice without over- 
lapping. Remarkably, most saddle point solu- 
tions obtained from more than one Z(2) vortex 
still exhibit well separated peaks after extremiza- 
tion. (We will call these saddle point solutions 
a" Z(2) saddle gas".) We can get various starting 
configurations by choosing different positions and 
orientations for the Z(2) vortices, and different 
U(l) subgroups. Starting from these configura- 



tions, we use deterministic cooling to get smooth 
configurations and then extremize further to get 
to the saddle point solutions. We used two dif- 
ferent extremization algorithms, one by Duncan 
and Mawhinney [3] and the other by Van der Sijs 
[5]. The final saddle point solutions are quite in- 
dependent of either the choice of extremization 
algorithm or the number of cooling steps initially 
taken. 

3. String tension 

Figure 3 shows the results of the measurement 
of the string tension of various SU(2) configura- 
tions. The string tension appears as a plateau in 
the Creutz ratio, given by 



C(R,T) = - In 



W(R,T)W(R-1,T -1) 
W(R - 1, T)W{R, T - 1) 



where W(R,T) is the average of the trace of a 
Wilson loop with size RxT. 

As shown in the graph, the string tension sus- 
tains its value for cooled Monte Carlo lattices, but 
was apparently broken for extremized lattices [3] . 
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Figure 3. The Creutz ratios as a function of R for 
a Monte Carlo calculation at (3 = 10.0 (o), cooled 
Monte Carlo generated lattices (□), extremized 
Monte Carlo generated lattices (O), and two dif- 
ferent 3-peak saddle-point lattices (A and *). 



For Z(2) saddle gas configurations, we see two 
kinds of behavior: some configurations exhibit a 
rising Creutz ratio, while others have a Creutz 
ratio falling off at large distances. The fact that 
most of the configurations exhibit either a ris- 
ing or falling of the Creutz ratio, not showing a 
plateau, shows that a plateau may be achieved by 
correct weighting. 

4. Eigenvalues 

The eigenvalues of harmonic fluctuations can 
be obtained by diagonalizing the matrix of the 
second derivatives of the action at the saddle 
point. 



(Ah 



d 2 S 
OAS A, 



where i, j are the indices for the gauge links. 

Because the size of the matrix we need to diag- 
onalize is large even for lattices of moderate size 
(3 x 3 x N 3 for N 3 lattice) , we should use some in- 
direct methods to calculate the eigenvalues. We 



use the Lanczos algorithm to tridiagonalize the 
matrix and use the bisection method to get the 
eigenvalues [6]. We have also implemented the 
spectral method [7] for comparison. For enough 
iterations, well separated eigenvalues agree be- 
tween the two algorithms. 

Table 1 shows the eigenvalues for a single and 
a double peak saddle. The first two columns are 
from a single peak saddle; the amount of extrem- 
ization differs. The second two are from a config- 
uration with two peaks. We didn't do any gauge 
fixing. 

Highly degenerate eigenvalues, especially zero 
modes from the gauge degrees of freedom, pose 
some difficulty getting all the eigenvalues using 
these methods. In the first column of Table 1, 
we can see many small eigenvalues of O(10~ 4 ). 
These are actually zero eigenvalues, mainly from 
the gauge degrees of freedom, spread due to the 
fact that the configurations are not extremized 
enough. In the second column, the eigenvalues of 
0(1O -4 ) have become £>(10- 10 ). This spreading 
happens to other eigenvalues as well, but is much 
more manifest for small eigenvalues. As you ex- 
tremize further, the amount of spreading becomes 
smaller and finally becomes less than roundoff er- 
ror, which is the case for the second column. 

Because of these close eigenvalues, we needed 
more than 6 x (size of the matrix) double preci- 
sion Lanczos iterations to get reliable results for 
the lowest eigenvalues (whether we could get all 
the eigenvalues for less extremized lattices using 
enough iterations is not clear.) For the spectral 
method, many eigenvalues in such a small interval 
make it impractical to resolve all the eigenvalues. 
Using the spectral method, we would need to it- 
erate for 



AA 

steps, where X max is the biggest eigenvalue, \ m in 
is the smallest eigenvalue and AA is the smallest 
difference between nearby eigenvalues. For the 
configuration in the first column, this is ~ 10 8 
iterations. 

As expected, these saddle point solutions have 
only a few negative eigenvalues, some of them 
degenerate, possibly from the symmetry of the 
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Table 1 

Lowest eigenvalues of Z(2) saddle gas configurations on a 12 3 lattice from the Lanczos algorithm. 



starting configuration. The differences between 
eigenvalues of configurations with 1 peak and 2 
peaks show that the quantum interaction with 
these Z(2) saddles can be important when we av- 
erage over configurations. 



5. Discussion 

Continuing study of the eigenvalue spectrum 
will show whether naive averaging of configura- 
tions will represent the partition function cor- 
rectly. More eigenvalue measurements and com- 
parison between the eigenvalues of the Z(2) sad- 
dle gas and Monte Carlo generated lattices are 
needed to study the validity of the starting con- 
figurations we have used. Also, measuring other 
physical quantities, such as the chiral condensate, 
on these configurations will give us an interesting 
perspective about the nature of confinement. 
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